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(D . Abstract 



We prove a useful identity valid for all ADE minimal S-matrices, that clarifies 
the transformation of the relative thermodynamic Bethe Ansatz (TBA) from 
its standard form into the universal one proposed by Al.B.Zamolodchikov. 
By considering the graph encoding of the system of functional equations for 
the exponentials of the pseudoenergies, we show that any such system having 
the same form as those for the ADE TBA's, can be encoded on A, D, E, A/Z2 
only. This includes, besides the known ADE diagonal scattering, the set of 
all SU (2) related magnonic TBA's. We explore this class sistematically and 
find some interesting new massive and massless RG flows. The generalization 
to classes related to higher rank algebras is briefly presented and an intriguing 
relation with level-rank duality is signalled. 
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1 Introduction 



In the recent years the understanding of the topological properties of the Renormal- 
ization Group (RG) space in two dimensions has undergone a remarkable progress. 
First of all, the discovery of the dissipative nature of the RG flows (the celebrated c- 
theorem has given a great insight into the problem. Then, for a large class of RG 
flows that show the property of integrability, the proposal by A.Zamolodchikov 
for the conjecture of a factorizable S-matrix corresponding to a certain perturba- 
tion of a Conformal Field Theory (CFT) by one of its relevant operators, allows 
a lot of non-perturbative information to be extracted. To give evidence that the 
conjectured S-matrix really describes the considered theory, one must extract from 
it information on the ultraviolet (UV) limit. This can be done using the proce- 
dure recently introduced in this context by Al.B.Zamolodchikov that goes under 
the name of Thermodynamic Bethe Ansatz (TBA) and whose original formulation 
traces back to Yang and Yang 0] . 

The TBA can be presented as a set of coupled non-linear integral equations 
driving the evolution of the Casimir energy of the theory on a cylinder along the 
RG flow exactly and non-perturbatively. In spite of their apparent complexity, they 
are often numerically integrable without using very heavy computer resources, for 
each point on the RG flow, and show the peculiar behaviour to be analytically 
solvable in the UV and infrared (IR) limits, thanks to transformations leading to 
sum rules of the Rogers dilogarithm function. 

The deduction of TBA equations directly from the S-matrix is easy only when 
the latter is a purely elastic (diagonal) one. In the most general case of non-diagonal 
S-matrix, one is lead to consider Higher level Bethe Ansatz to get a TBA system out 
of it. This is a formidable task in many cases, and an alternative strategy would be 
welcome. In a beautiful and stimulating piece of work Al.Zamolodchikov noticed 
that the TBA system for purely elastic scattering matrices related to A,D,E Lie 
algebras (i.e. those previously treated in P), can be suitably transformed in a form 
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where it appears a clear encoding on the Dynkin diagram of the related A,D,E 
algebra. 

Reversing the strategy, one can draw a diagram, set up a TBA on it, and 
compute formally the central charge of the CFT (with action, say, Scft) describing 
the UV limit, as well as the conformal dimension of the perturbing operator $. It 
is then reasonable to conjecture that the theory described by the action 



S = Scft + \j d^x^{x) (1.1) 



put on a cylinder, has a Casimir energy driven along the RG flow by the proposed 
TBA. At this point other information can be extracted to test further this conjecture 
and to explore, for example, the IR limit of such a theory. Following this very 
productive attitude, in some recent papers [0, ||, |^, |10|, [ill a lot of RG flows have 



been studied, and even some new discovered [12 



While in the diagonal TBA, energy terms in the equations, one for each particle 
type, were naturally attached to a// the nodes of the corresponding Dynkin diagram, 
in the non-diagonal case the structure of the encoding is a bit different, as there 
are a quantity of nodes with no energy term attached. These nodes correspond to 
what in Bethe Ansatz literature are called magnons, i.e. fictitious particles with no 
mass and no energy, whose unique task is to exchange internal degrees of freedom 
(colors) between the physical particles of the theory. We shall often refer to such 
TBA systems as magnonic. 

In the present paper we attach the program of systematically exploring this 



large class of TBA systems. In |TT] it was realized that for a large set of models, the 
encoding of TBA is natural on a certain kind of "product" of two Dynkin diagrams, 
one for the physical particle structure, the other for the magnonic structure. In 
particular, when there is only one mass in the spectrum (so that the "physical" 
Dynkin diagram is Ai), the magnonic TBA are encoded on a single Dynkin diagram. 
To be more precise, in |jll| only the case of An magnonic diagrams was explored. 
One can ask if the TBA encoded on, say, D„ or diagrams have any meaning. 
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and even further, if there are some other class of graphs, not Dynkin diagrams, that 
can encode the magnonic structure of TBA. In the present paper we give an answer 
to these questions for the case where the "physical" Dynkin diagram is Ai. It turns 
out that the magnonic TBA has some basic properties (the so called Y-system, see 
below) in common with the set of diagonal TBA's encoded on ADE diagrams. 

For the latter, this Y-system is generated thanks to a useful identity on S- 
matrices that we prove in sect. 2 (after review of the needed formalism of [p!3| ). 
Then in sect. 3 we use this identity to give a proof of the transformations that lead 
Al. Zamolodchikov to his universal form of TBA |^. We also get the Y-system (i.e. 
a system of functional equations to be satisfied by the solutions of TBA). Sect. 4 is 
devoted to the proof that a Y-system of the kind found in sect. 3, can be encoded 
only on ADE Dynkin diagrams (plus the tadpole diagram T„ corresponding to 
a folding of A2n)- We show that a class of magnonic TBA's generalizing those 
proposed in 0, H, ^, |l^ has a Y-system that simply maps into the previous one, 
thus allowing to extend the ADET classification to this case. Next, in sect. 5, we 
explore systematically the whole set of TBA's thus proposed, trying to identify their 
UV limit, the perturbing operator, and, when applicable, the non-trivial IR limit. 
While many of the flows thus described were already known, some new appear, 
especially in the study of the -^6,7,8 and T„ cases. We end in sect. 6 by commenting 
on the generalizations when the "physical" diagram is not Ai, on the possibility to 
envisage a general scheme for all TBA's and putting a remark on a still mysterious 
relation with level-rank duality in CFT. 



2 ADE S-matrices: a useful identity 

We briefly summarize some basic facts about purely elastic scattering theories which 
will be useful in the following. For our purposes it is convenient to start from Dorey's 
approach to the ADE S-matrices |TB[. A (1 + 1) dimensional purely elastic scattering 
theory has a factorizable and diagonal S-matrix. Factorizability means that the 
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scattering amplitudes of any number of particles can be written as products of two- 
particle amplitudes. Therefore, the scattering of particles a and b is described by 
the two-particle scattering amplitude Sab, which is a function of the relative rapidity 
^ab =\ (^a — (^b \- A simple pole of Sab at 9ab = iU^b ^^e direct channel indicates 
that there exists a bound state c of a and b whose mass is: 

'ml='ml + ml + 2 rua rrtb cos(f/^J (2.1) 

and the scattering amplitudes with any particle d must satisfy the bootstrap equa- 
tion: 

ScM = Sad{e + iul)SM{e - lU^,) (2.2) 

where tJ^b = n — [7^^. For the conserved charges, the bootstrap equation leads to 

e-isul,qa ^ e^sus,qb ^ (2.3) 

In the so called ADE scattering theories the fusing angles U are all integer mul- 
tiples of ^, /i being the Coxeter number of the Q = ADE Lie algebra (of rank r) 
associated to the theory. It turns out that nontrivial solutions to the conserved 
charge bootstrap only occur if the spin s modulo h, is equal to an exponent of Q. 
Furthermore, each of the r particles in the theory may be assigned to a node on 
the Dynkin diagram of Q, in such a way that the set of conserved charges of spin 
s, when assembled into a vector = (q'"S Q'"^? ••Q'"''); forms an eigenvector of the 
incidence matrix G of the Dynkin diagram of Q with eigenvalue 2cos(^). Thus, 

Gq, = Xsqs , A, = 2cos(e,) , 0, = — (2.4) 

n 



Notice that for s = 1 eq. (p.4|) gives the masses of the particles in the theory, thus 
showing that they are organized in the so called Perron- Frobenius eigenvector of G, 
namely in the eigenvector ipc corresponding to its highest eigenvalue. For a matrix 
with non-negative integer entries like G, ipc turns out to be always unique and has 
all non-negative components. 
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Let $ be the root system of Q. Its Weyl group, i.e. the group of all reflections 

i CX Xi 

Wa{x) = X - 2 ) ' . a , (2.5) 
[a, a) 

that map $ into itself, is generated by the subset of the reflections associated with 
the set of simple roots 11 = {ai, ...0;^}. The Coxeter elements of the Weyl group 
are elements of the form WaiWa2 ■ ■ - War- Sphtting 11 into two subsets of orthogonal 
roots: 

U = (ai, a2 ■ ■ ■ ak) U {Pi, p2 ■ ■ ■ Pr-k) , {a^, aj) = {pi, pj) = 2Sij (2.6) 
one defines 

W = Wa^Wa2 ■ ■ ■ Waf,Wf3^Wf32 ■ ■ ■ Wp^_^ (2.7) 

w is called a Coxeter element and has period h. The group generated by w, is 
therefore isomorphic to Zh- Let {a, (3} be the dual basis to the simple root {a,P}, 
\s = 2cos{9s) an eigenvalue of the incidence matrix G not equal to 0, and Qs the 
corresponding eigenvector. Defining 

as = J:C<^^ bs = J2C^ (2.8) 



following the arguments of fl^, one can see that, for the simple roots 

{at, as) = qT {Pi. a^) = (a„ 6,) = (A, 6,) = gf' (2.9) 

and we can define a projector Pg into the two-dimensional subspace spanned by 
and hs 

Ps{a^) = qTas Ps{-Pj) = q^K (2.10) 

where {cis,—hs] are dual to {as,hs] in that subspace. Since and have equal 
magnitude, this implies that the projections of the simple roots have lengths pro- 
portional to the components q\ of the eigenvectors of the incidence matrix. The 
Coxeter element w acts in each subspace as a rotation by 29s. Hence we have 
(introducing a complex notation in each invariant subspace) 

Ps{w^a,) = g^"'e^(2p+i)e. (2.11) 



Ps{wP{-pj)) = -gf^e^(2p)e, ^2.12) 

In this formulation the general expression for the ^af,— matrix element in the ADE 
scattering theories are: 

(a) Particles a and b of type a 

^a6= (2.13) 

(b) Particles a of type a , 6 of type /3 

Sab = n {2^}?'"''''"'^ (2.14) 

p=0 

(c) Particles a and b of type jS 

^a.= nW-l}?"'"'"'^ (2.15) 
p=0 

where 

{a;}+ = (x- l)+(x + 1)+ (2.16) 
, , I , i7ix\ 

We use this formalism to prove a useful identity: 

Sab {O + ^) Sat (^^ - ^) = n Sacief"^ , M (2.18) 



The proof goes as follows. By using eq. ( |2.11| , p.l^ ) it is possible to rewrite equation 

(U) 

Gabql = e-q: + e—q: (2.19) 



as 



aa + UJ ^aa= -^Gabl3h , -uPb- Pb = J2GbaO:a (2-20) 

b a 



Then we have for case (a): 

Sab {9 + ^) 5., {d-'^) = n {2P + 1}?-"''"^^ n {2fc - 1}?-'"'"^) (2.21) 



p=0 fc=0 
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or, after a rescaling k ^ p + 1 (and using the property {2h + x}+ = {a;}. 



Sab{d + -j^^ Sa 



and using the identity ( |2.2q ) 



- =n{2p+i} 

^ p=0 



(aa,w'P(-/36-ii;/36)) 



(2.22) 



ITT 



Sab ( ^ "/^ ' '^'^^ 



^-x) = n{2p+i} 



h-l 



J2k Gbk{oia,w''ak) 



USacief"^ (2.23) 



p=0 



The other cases (a) and (b) ( p^.lSj ) go through the same way, and so ( |2.13[ - PT5| ) does 
indeed provide a set of functions which obey the S-matrix eq ( p.l8| ). 

For 6 = identity (|2.18|) must be carefully treated, as it becomes apparent by 
taking its logarithmic form 

log Sab + X j + log Sab [d-j)=T. Gbc log SaM - 2t7rQ{d)Gab (2.24) 
The term proportional to the step function 

if X < 


1 if X > 

has to be introduced to take into account the correct prescription for the multivalued 
function logx. According to (|2.24|), formula (|2.18D must be corrected as follows (if 



Q(x) = lim 



1 1 X 
— I — arctan — 

2 vr e 



I if X 



(2.25) 



we want to include the point 9 = 0) 

Sab + -j^^ Sab (d 



^) =n^-(^)^^=e-2^-^-®(^) 



(2.26) 



The corrective exponential term is 1 for all values 6' 7^ 0, while at = corrects the 
r.h.s. of the identity to be compatible with the fact that the l.h.s. becomes for 9 = 
the unitarity constraint on the matrix 5, while the values «S'(0) must reproduce the 
correct statistics of the system. We shall appreciate the deepness of this corrective 
term in next section, where we relate it to the TBA equations. 

We would like to emphasize that (p.26|) often gives relations equivalent to some 
of the bootstrap equations ( |2.2| ). The deep interrelation between our identity and 
the bootstrap certainly needs more investigation. 



3 Universal form of TBA 



In it has been proposed, to recover the information on the ultraviolet (UV) 
limit of the theory defined by the matrix Sab, to use the so called Thermodynamic 
Bethe Ansatz (TBA), which is a set of non-linear coupled integral equations driving 
exactly the Casimir energy of the system (on a cylinder of circumference R) along its 
Renormalization Group (RG) flow, thus allowing the determination of the effective 
central charge c = c — 24Ao of the UV theory, where the lowest conformal dimension 
Aq is for unitary theories (for which then c = c), and negative for non- unitary 
theories. Putting Ua = rriaR cosh 6 (the so called energy term), the TBA system 
is a set of equations in the unknowns Ea (often called pseudoenergies) , having the 
general form 

iya{0) = Baie) + ^ Y.i<f>ab * log(l + e-^))(0) (3.1) 

where 0a6 = ""^^^ log ^ab and the * stands for the rapidity convolution 

/ + 00 
d9'A{9 - 0')B{d') (3.2) 
-oo 

From the solutions to this system, the evolution of the vacuum energy E{R) = 
— ^g^^^ along the RG flow can be followed by use of the equation 

3 /- + 00 , ^ 

ciR) = - E / ^a(^) log(l + e~''^^'^)de (3.3) 

which, in the i? — > limit, turns out to be expressible, after some manipulation 
involving the derivative of eq. p.l| ) (see for example ref. Q), in terms of Rogers 
DilogarithmQ sum rules 

5 = 5(0) = 4 (7^1 (3-4) 
with ya given by the solutions to the algebraic trascendental equation 



/■+00 (jQ 

y„ = n(i + iM)^''^ , Nab = - -Me) (3.5) 

^ J —00 ATT 



00 



"'^see for example |l4| for a definition and for properties. 



that can be deduced from eg. ( |3.1| ) in the hmit i? — > 0. 

In a nice recent piece of work P| ALB.Zamolodchikov proposed a transformation 
of TBA equations for ADE diagonal scattering showing in a clear form their relation 
of the set of integral equations to the ADE Dynkin diagrams. This transformation, 
leading to what is now known as universal form of TBA, is based on a remarkable 
matrix identity quoted in |^ 



1 T 



'■>ab 



G. 



ab 



(3.6) 



2 cosh(7rA;//i) 

where (pabik) stands for the Fourier transform of ipabid) 

<Pab{k) = / dd<f)ab{e)e"'' (3.7) 

J —CO 

As in there is no explicit proof of this identity, we give here a proof based on 
our identity ( p.24| ). We shall also explain in some more detail how the universal 
form of TBA can be deduced out of it, and also how one can obtain a system of 
functional equations which is also given in [^] and is very useful in order to extract 
further information on the TBA system. 

First of all, let us derive Zamolodchikov's identity from ( |2.24|) . Take the deriva- 



tive of eq.( |2.24| ) and define as usual (pabifi) = —i-^^ogSabid) 



J2Gaac{0)-27lS{9)Gab 



(3.8) 



and Fourier transform this equation {k is the momentum corresponding to 

A;7r\ 



or 



2 cos 



<Pab{k) 



h 



^ab 



{k)=Y.G,,cPac{k)-2TxG, 



ab 



-2tx 



G 



(3.9) 



(3.10) 



ab 



This equation is trivially equivalent to the matrix identity ( p. 71) , but form (p.9|) is 
even more useful for our purposes. Notice that (|3.10|) computed in A; = recovers 



a well known identity P, [15 



Ar = A(2-A)-^ 



(3.11) 



and this helps to transform eq.( p.5| ) into the more appeahng form 

yl = Ui^ + ytf^' (3.12) 

6 

Fourier transform eq. (|3.1|) , then muhiply both sides by Sab — R{k)Gab, where 
-^(^) ~ 2cosh(7rfc/fe) ' finally use ( p.9| ) to recast the TBA system in the following 
universal form 

M^) = Baie) + ^ E GablVh * i^b - log{l + e'^ym (3.13) 

b 

In this form, the TBA is explicitly fixed once the diagram whose incidence matrix is 
G is given. The universal kernel (fh, which is (up to 27r) the Fourier antitransform 
of R{k), depends only on the Coxeter number h of G 

= jJkw) 

Notice that in the R ^ limit eq.( p.l2| ) is directly obtained instead of ( p.5| ). 

Now let us consider eq. (|3.1| ) for 9^9 — in/h and for 9^9 + in/h. Summing 
up and subtracting eq. (|3.1|) calculated in 9 and multiplied by G and using ( p?9|) we 
get 

Uai9+i7r/h) + Uai9-i7r/h)-Gabi^bi9) = 6ai9+in/h) +ea{9-m /h)-Gabeb{9)-GabLb 

(3.15) 

finally using the identity (^^ ) for s = 1 

Ya (e - ^) Ya {e + '^)= n(i + Y,{9)r^^ (3.16) 

Where Ya = e^"'-^-*. This system (that we call Y-system in the following) is ex- 
tremely important, as commented by many authors, as it seems to encode even 
more information on the system than the usual TBA. 

First of all notice that the stationary solutions of this system (i.e. those who 
do not depend on 9) are exactly the Ua appearing in eq. ( |3.12| ), which are the basic 
tools to extract the UV central charge. 
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An '■ 



En 



1 2 i n-l n 



1 2 3 n-2 



1 2 3 n-2 n-1 



a = n + 1 — a, a = l,...,n 

for n even a = a, a = 1, . . . ,n 

{a = a, a = l,...,n — 2 
n = n — 1 

for n = 6 : 1 = 5, 2 = 4, 3 = 3, 6 = 6 
for n = 7, 8 : a = a, a = 1, . . . , n 

n < 8 



T • 



1 



3 n-1 n 



^2n/^2), a = a, a = l,...,n 



Figure 1: A,D,E,T diagrams: the numbers show the labelling of the differ- 
ent nodes. On the right the particle-antiparticle relations between nodes are 
shown. 



Moreover, as stressed in the Y-systems encoded on Dynkin diagrams show 
a remarkable periodicity 

Y,{e + Pzn) = Y^{e) , P = ^ (3.17) 

h 

where a represents the antiparticle of a (see fig.l). 

This can be shown (along the lines of 0) to be in relation with the conformal 
dimension of the perturbing field, via the formula 

A = 1 - for E„ , A = 1 - ^ for T„ (3.18) 

This allows to extract in a simple way the parameter A, characterizing, together 
with the central charge c, the action of the theory. The explicit proof of this 
periodicity relies on successive substitutions inside the Y-system of the functions 
Ya computed at different points. This is in general a very cumbersome task, even 
for the most simple cases. We do not know of a general proof of the periodicity, 
better we used a simple computer program to test it up to 16 digit precision for all 
Dynkin diagrams up to rank 50. 
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4 Y-systems on general graphs 

In the previous section we got the Y-system from ADE massive scattering theory. 
We can ask if such a system can be more generally defined on graphs other than 
the Dynkin diagrams, thus allowing generalizations of the ADE scattering theories. 
Here we prove that any Y-system of the form l \3. 1 (\ ) where Gab is a matrix with non- 



negative integer entries, can allow for stationary non-negative solutions (Va ^ ^) 
only if Gab is the incidence matrix of an An, Dn, -£^6,7,8 Dynkin diagram or the 
"tadpole" graph Tn = v42„/^2- This set of graphs, shown in fig.l, defines the set 
of all square matrices with non-negative integer entries whose norm (the highest 
eigenvalue, corresponding to the Perron- Frobenius eigenvector) is strictly less than 



2 [|T6|. The requirement of existence of a stationary solution is a must to have a 
well defined central charge for a system described by a TEA of form ( p.l3|) . 

As Ea are real functions, e"^"^ = Ya > 0. In particular, for solutions not depending 
on 9, ya > 0. More precisely, it can never happen that ya = 0, otherwise in eq.( |3.12|) , 
at least one of the factors on the r.h.s. must be zero, implying one of the yb to be 
— 1, in contradiction with its positivity. Therefore we assume ya > for all a. 

Defining Xa = \/ya + 1, eq.( p.l2|) becomes^] 



xl-\ = YixX-' (4.1) 

b 



with > 1. Moreover, we can pose Za = logXa (and then Za > 0). This allows to 
write the logarithm of the l.h.s. of (|4.1|) as 

1 ~ 



log(x^ - 1) = 22;, + log ( 1 - ) < 2za (4.2) 



Therefore, from (WA 



^ GabZb < 2Za (4.3) 

b 

Za can be decomposed in the base of eigenvectors of G and, having all positive 
components, it has a positive and non-zero projection on the Perron-Frobenius 



^Thc following argument has been suggested to us by F.Gliozzi 
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eigenvector. A projection of formula ( [4 .31 ) on the Perron-Frobenius direction simply 
shows that Ap < 2, where Ap is the Perron-Frobenius eigenvalue. This bound is 
known to select the incidence matrices of the graphs A, D, E, T drawn in fig.l. 
Corresponding to the Y-systems of the form ( p. 16 ) we have seen that there are 



TBA systems of the (universal) form ( p. 131) . These are nothing but the whole set 



of TBA's studied in the paper of Klassen and Melzer (including the case they 
call that corresponds, in our notation, to the T„ diagrams). We refer to 
for a complete description of the identification of the models at UV, and of their 
perturbing operators. 

We notice, however, that the proof of classification of Y-systems we have given 
is absolutely independent of h. Other choices of the parameter h, where it no more 
plays the role of Coxeter number, can, in principle, lead to sensible TBA systems. 
One such choice, on which we shall concentrate in the following, is the magnonic 
TBA proposed by ALB. Zamolodchikov [0, H, § to describe RG flows of minimal 
models perturbed by their least relevant operator 0i3. This TBA has the general 
diagrammatic form 

2tt ^ cosh d 



The rationale under this form of TBA will appear later |]TT|. Here the terms Va are 



zero on all nodes but one or two (labelled /c, / in the following formula) 
, S^mR cosh 9 for massive case, A < 

iyaie) = { (4.5) 

= + 5^6 ^) for massless case, A > 

Taking this equation for 6 summing it to the same equation for 6 ^ 6+^ 

and taking into account the pole structure of the kernel 0=1/ cosh 6 and the fact 
that if Ua are given by eq.([4.5|) then z/a(6' + in/2) + Uai^O — 111/2) = 0, one arrives at 
an equation having the same form of ( p.l6|) , where now Y^{9) = e-^"(^) (notice the 
different sign in the exponential). The stationary solutions to this system again are 
the basic tool to compute the UV central charge. 
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The surprising fact is that in this case the hnes of reasoning that lead to the 
A, D, E, T classification of Y-system (|3.16| ) apply as well, therefore such a clas- 



sification holds for this magnonic structures too. Next section is devoted to the 
exploration of all possible such magnonic systems on A, D, E, T. 

5 ADET magnonic TBA 

To begin this section, we would like to emphasize some general rules on the dia- 
grammatic approach to the magnonic TBA ( |4.4| ), that we use in the following. Then 
we briefly comment on non-perturbative terms, and finally present our results on 
the systematic exploration of all ADET cases. 

5.1 cuv massive models and dilogarithm sum rules 

For massive models, imagine that the single massive term Ua different from zero is 
put on node k. Then to compute the UV central charge consider the full diagram 
G and the diagram G' = G — {k} where the node k and the links emanating from 
it have been deleted. Then, following the arguments in p, ||, |^, |, ^, we have the 
following general rule 

Cuv 



ciR instead is zero, as the model is massive and has a trivial IR point. 

Formula ( |5.1| ) is easily computed by resorting to the dilogarithm sum rules 
quoted in and remembering the following identity on Rogers dilogarithm 

\l + x) 6 \l + x) ^ ' 

5.2 cuvt ciR and parity issues for massless models 

For massless models, call R the node on which the left mover z/ = is put, and 
L the node on which the right mover u = '^e~^ lies. Then consider the diagrams 
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G' = G — {L} and G" = G — {L, R}. The UV central charge is given by 

cuv = ^\y:-y:\L\-^\ (5.3) 



7r2 



and therefore coincides with the calculation for the corresponding massive case, 
where the mass is put on the same node as the left mover. The IR central charge 
instead is given by 

Notice that the position of the right and left movers, i.e. the choice of nodes 
L and R is not arbitrary. Parity invariance of the vacuum requires that the TBA 
system must be invariant under exchange 9^—9 and therefore R and L must be 
interchangeable on the diagram without changing the TBA structure, i.e. they must 
lie on nodes symmetric with respect to some Z2 symmetry of the diagram. This 
is a strong constraint on the possible massless flows. For example Ey^^^Tn do not 
possess any Z2 symmetry and it is not possible to write a sensible TBA describing 
a massless flow on them. This is in connection with the parity of the perturbing 
operator. If the operator is even, two different behaviours are in general expected 
for different signs of the perturbing parameter, one being massive and the other a 
massless crossover to a non-trivial IR theory. Conversely, for parity odd operators 
the sign of the perturbing parameter can always be readsorbed in the operator and 
does not affect the (massive) behaviour of the perturbed theory. 

5.3 Periodicity of Y-system and conformal dimension of the 
perturbing operator 

The conformal dimension A of the perturbing operator can be deduced from the 
periodicity of the Y-system and is independent on the choice of the particular nodes 
where masses or left-right movers are put. Of course, as the role of h is changed, 
the periodicity also gets some modification. Eg. ( p.lTp still holds, but now 



P = ^ (5.5) 
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Moreover, the symmetry on the diagrams for a ^ a is now destroyed by the asym- 
metric choice for Ua. When a 7^ a in fig.l, the real periodicity is doubled. One 
should be careful, however, that this prediction can be affected by some selection 
rule on correlation functions of the perturbing operator at criticality coming from 
symmetries of the conformal fusion rules governing the UV CFT. A careful analysis 
of all the cases leads anyway to a formula like ( p.l8| ). For reader's convenience, we 
list here the results for all diagrams 

A„,: A = l-^ E,: A = l 

I^n: A = l-i Er: A = ^^ (5.6) 

T„ : A = 1 - 1^ Eg : A = § 

5.4 Non-perturbative terms in the Casimir energy 

Once the UV and IR behaviours are identified, the expansion of the Casimir energy 
in terms of R is an issue. This contains both perturbative (in g = R^""^^) and non 
perturbative terms. The non-perturbative contributions can be computed along 
the lines of 0, |^, 0, |^, 0, and amount or to a bulk term proportional to R^ or to a 
logarithmic term proportional to R^ log R. This latter appears when the incidence 
matrix G is not invertible. In other words, the scale function F[R) = for 
Eq,E8, has the general form 

12 2n 

while in the other cases Dn, A2n+i, Eq the logarithmic bulk term appears 

F{R) = -- + L,\og{mR)iy-^^ +Y.Ug'' (5.8) 

where the coefficient pertains only the node k where the mass (or the left mover) 
is put, and can be elegantly expressed by considering the eigenvector q correspond- 
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ing to the null eigenvalue, i.e. 



An odd- g=(l,0, 
Dn odd- q = (0, 0, 



i,o,i,---) 

-,0,-1,1) 



D. 



E- 



7 • 



n even 



g= (1,0, 0,0, 0,-1,-1) 

g« = (0,0,---,-l,l) 

g(2) = (...,0,2,0,-2,0,1,1) 



(5.9) 



and correspondingly 



A. 



■n odd - 



T — ^ 




E7 



Lk 



(5.10) 




ilk? 



D, 



n even 



Lk 




5.5 Model identifications 

In the following we summarize and comment the results concerning our exploration 
of all possible ADET magnonic TBA structures. We divide the list according to 
the Dynkin diagram, and for each Dynkin diagram we first put a massive energy 
term Vk = mr cosh.6, and let k vary along the diagram up to exhaustion. Then, 
if the diagram presents some Z2 symmetry k ^ k, we put z/^ = mRe^ /2 and 
vj. = mRe"^ /2 and let k vary on all nodes of the diagram with non trivial image 
under this Z2. 

5.5.1 An case 

The work concerning An Dynkin diagrams has already been done completely by 
Al.B.Zamolodchikov in the series of works 0, ^, 0. For reader's convenience, we 
summarize here his results. To take advantage of the Z2 symmetry of the An 
diagram, it is convenient to put n = k + I — 1 and consider A^+i^i. The z/^'s are 
chosen as in ( [4.5| ). The resulting flows start at the (fc, /)-th SU{2) coset model, 
perturbed by its operator of dimension A = 1 — -j^^-^ and when A < evolve 
into the massive theory described by the non-diagonal S-matrices of Ahn, Bernard, 
LeClair |]T^, while for A > the flow is massless and at the IR limit reaches the 



{k - I, l)-th 5f/(2)-coset. 
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5.5.2 Dn case, mass on the tail 

Putting a mass term i/^ = mR cosh 6 on tlie k-th node of tlie tail {k = 1, ...,n — 2) 
of a Dn diagram, one can describe a massive RG flow whose UV limit has central 
charge 

(5.11) 



k + 2 

and therefore lies on the k-th critical line of those described in This crit- 



ical line, in turn, can be seen as the UV limit of a fractional super-Sine-Gordon 
theory , and the perturbing operator can therefore be identified with where 

$ = ,/,^ : eTfe"^ : (5.12) 

ipi being the generating parafermion of dimension 1 — 1/ k and (j) a free massless 
bosonic field, so that the vertex operator in ( p.l2[ ) has dimension /S'^/Stt. The 



dimension of $ must fit the value predicted by the periodicity of the Y-system 
( ^.6|) . Notice that the value of c depends only on k (the node where we put the 
mass) and not on n, the rank of the diagram. Hence, for each k there are a sequence 
of points on the critical line labelled by n. The identification is done by comparing 
the dimension of the perturbing operator as predicted by the periodicity of the 
Y-system, namely 1 — 1/n and the dimension of $ as described above. This yields 

(3' 1 I . . 

— = --- 5.13 

8tt k n ^ ^ 



Notice that this result for k = 1 was known to Al.Zamolodchikov, as quoted in |[T2 
This allows to identify the S-matrix of the perturbed massive theory as 

S = Sk®SsG{l-l^) (5.14) 



■^By considering the ground state TBA we can not distinguish a model from its orbifold (this 
could be done by considering excited states TBA, instead), therefore here and in the following we 
"identify" models differing one from the other by some orbifolding procedure. For example, we 
speak of a single critical line at c = 1, while it is known that a more subtle analysis shows [p^ 
two lines which are one the orbifold of the other. 
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where Sk stands for the Bernard LeClair I^O] k-th minimal model + 0i3 S-matrix 
and SsciP"^ /^T'') means the Sine-Gordon S-matrix at coupling p. 

Notice that our TBA for such S-matrix is in agreement with the recent observa- 
tion in about the gluing (at the "massive" node) of diagrams pertaining each 
factor in a TBA corresponding to a tensor product S matrix. 

Finally we would like to signal that when k = n — 2, i.e. when the mass goes on 
the bifurcation point of the diagram, we have the N = 2 supersymmetric point of 
the corresponding critical line, as stressed recently by Fendley and Intriligator 



5.5.3 Dn case, mass on the fork 

The other possibility is to put the mass on the node n (or equivalently n — 1). This 
case has been analyzed by Fateev and Al.Zamolodchikov |[T^. We report it here for 
completeness. The UV central charge calculation gives 

^uv = ^ (5.15) 

n + 2 

thus showing a dependence on n in this case. This turns out to be the central 
charge of the celebrated Z„, parafermionic models, SU{2)n/U{l) as coset models. 
The perturbation, as usual, is identified by the periodicity of the Y-system to have 
A = 1 — ;i, and it is therefore identified with the operator tpi{z)ilji{z) + ?/'|(z)?/'|(z), 
where ipi is the generating parafermion. This perturbation is parity even, hence we 
expect it to be sensitive to the sign of the coupling A. 

Indeed the Dn Dynkin diagram has an evident Z2 symmetry exchanging n with 
n — 1 and allowing the definition of a massless TBA flowing in the direction opposite 
to the previous one. In this case put a left mover on n and a right mover on n — 1. 
The UV central charge is the same as before, but the IR one is now given by 

6 

^iR = '^- 7 — rrvt — (5-16) 

{n + l){n + 2) 

thus giving evidence of a deeply non-perturbative flow between the Z„-parafermion 
model and the n + 1-th minimal model. This result shows how the approach of 
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conjecturing a TBA and then trying to identify the UV and IR hmits by use of the 
Y-system is a very effective one: it can give evidence of highly non-trivial and even 
unexpected results on the structure of the RG space in two dimensions. 



5.5.4 En case: mass on node 1 

This is perhaps the most intriguing case of our analysis. The three cases yield the 
following values of the UV central charge and of the perturbing operator conformal 
dimension 

E, c=l A = f 

E7 c=g A = A (5.17) 
Es c=l A = i| 

To identify the sequence of models giving the UV limit, it is interesting to complete 
this table by extending the En diagram to < 6 by taking E^ = D^, E4 = A^, 
E^ = A2 (B Ai. The second Ai factor in the last case is a pure magnon decoupled 
from the theory and it drops. Only the first factor is relevant. This allows to extend 
the previous table with the additional cases 

Es A = | 

c=f A = f (5.18) 
E, c=l A = | 

What is peculiar with all these UV models is that they share the property to be 
invariant under a generalized parafermionic algebra with grading (i.e. fusion 
rules ipi X ipj = ipi+j mod k) and generating parafermion ipi of dimension 1 + 1/K. 
These algebras have been called SZ^ in . Here K is related to n oi Enhy K = 



n—1. It is expected that such theories are the UV limit of the S-matrices proposed 
in [^] having Z;^-exotic super symmetry. The surprising fact is that, as shown 
in p4| (see also [^), with the hypothesis that no operator of dimension 1 appears 



as secondary of the identity, this series of algebras truncates at < 6. To be more 
precise, the SZ2 algebra is the = 1 superconformal algebra, generated by a field 



of spin 3/2, then SZ3 is the spin 4/3 algebra of Fateev and Zamolodchikov ||27 
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the SZ4 algebra describes a model on the c = 1 critical line invariant under a 
symmetry generated by a spin 5/4 field, the two remaining cases are degenerate: 
SZ5 = Z5 and SZq = ^2(8)^3, where Zn are the usual parafermions. Notice that the 
central charges of the models described by these algebras (for K = 2,3 the bottom 
models of the relative series of minimal models) are exactly those arising in our 
TBA computation. It is then tempting to identify the off-critical versions of these 
models (perturbed by the operators indicated above) with the scattering theories 
having the Bernard Pasquier S-matrices. However a first question immediately 
arises: what corresponds as UV to S-matrices with K > 67 This is up to now 
an open question. Secondly, we could for a while rejoice seeing that the truncated 
series of SZ^ corresponds to an En series, which is truncated too. However the two 
truncation are not at the same K. The case of seems not to enter this framework. 
The value of 3/2 for the central charge for Eg suggests that we are dealing with a 
model at a specific point of the super-Sine-Gordon line. If we try to identify this 
model by use of the perturbing operator dimension (as before for the -D„'s) we get a 
theory with p^/Sn = 7/16 that can be shown to possess a parafermionic symmetry 
generated by a ipi of dimension 8/7! (now there is an operator of dimension 1 as a 
secondary of the identity, hence the hypothesis of [|^, ^ do not apply here). What 
are the generalizations of this series for, say, parafermions of dimension 9/8, 10/9, 
etc.? This is an intriguing problem, and although its solution lies out of the scope 
of the present paper, we intend to return on this point in future. 

To conclude, we notice that the Ej^g diagrams do not possess any Z2 symmetry 
and it is not possible to define massless fiows on them. However, for the Eq diagram, 
we can transform the mass on the first node to a left mover and put a right mover 
on the symmetric node 1 = 5. This shows that the perturbation of the Z^ model 
by its second energy operator 62 of dimension 6/7 is even. For negative values of 
the perturbing parameter it fiows to a massive scattering theory described by the 
aforementioned Bernard- Pasquier S-matrix, while for A > this defines a massless 
theory fiowing to an IR limit that can be easily be computed to have c = 1. 
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Figure 2: Eqj^s Dynkin diagrams: the numbers are 
the values of the UV central charge when the massive 
energy term is put on the corresponding node. 



5.5.5 En case: mass on other nodes 

For En and mass terms on nodes other than 1, the results of the calculations of 

cuv have been encoded, for reader's convenience, in fig.2. We do not enter in much 

detail on the identification of UV models for these cases. Most of them can be 

identified with tensor products of mutually non-interacting minimal models. 

We discuss a single case which can have some interest by itself. Putting a mass 

term on node 2 of the Eq diagram we get cuv = ff , which corresponds to the tensor 

product of two copies of the m = 7 minimal modelQ The A = | perturbing field is 

realized as the tensor product of operators of dimensions 3/28 and 3/4 respectively. 

The perturbation happens to be parity even as one can figure out from the known 

parity of operators in m = 7 minimal model. This is in agreement with the Z2 

symmetry of the Eq diagram. Therefore, by replacing the mass on node 2 by a left 

mover and putting a right mover on node 4 we recover a massless flow between the 

aforementioned model and an IR limit with cjr = 81/70 that we can identify with 
"^The other possible identification with the Z26 parafermion is ruled out by the non-existence 
of an operator of dimension 6/7 in its Kac-table 
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the m = 5 model of the = 1 superconformal series. 



5.5.6 Tn case: a new series of non-unitary massive flows 

Another quite unexpected result concerns the tadpole diagrams T„. Here it is 
convenient to introduce the parameter p = 2n + 1 = 3, 5, 7, .... Notice first of all 
that the T„ diagram has no Z2 symmetry, so we can expect pure massive flows only, 
and the perturbing operator will be odd. Put a mass term Ui = mR cosh 6 on node 
/ of T„. The central charge computation gives 



l + 2\ p{p + 2l) 

The central charge computed here is, to be precise, an effective one c = c — 24 Aq, 
where Aq is the lowest conformal dimension in the UV model (negative if the model 
is non-unitary). Taking into account this fact, it is possible to identify the UV 
models with non-unitary '^^^^'(f)'^^^^^^' cosets, with / integer and k = | — 2. The 
perturbing operator with A = turns out to be the usual 0i,i,3 field. The first 
series / = 1 is given by the non-unitary minimal models perturbed by their 

013 operator. The second series is supersymmetric, the models are One 
of these models also belongs to the minimal series, namely the first SM^j = M7 12. 
This gives therefore also a description of the flow of the theory Mj^u + ^4,5. 

The S-matrices of these models are a tensor product of minimal Bernard LeCLair 



S-matrices times Smirnov reductions of the Sine-Gordon S-matrix for fractional 
value of k 

S = Sk(^Si , leZ , keZ + ^ (5.20) 
This also is in agreement with the "gluing" procedure suggested in 



6 Conclusions, generalizations and final remarks 

We have explored the whole class of magnonic TBA's whose Y-system is of the form 
( |3.16 ). Of course, this is far from being the most general case. In [0 it has been 
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shown how the TBA's for higher coset models perturbed by 4>id,id,adj organize in a 
nice way in terms of two Dynkin diagrams, one pertaining the physical particles 
(call it G) and one the magnons (call it H). The general TBA for coset models has 
the universal form 

+ ^0. * (e Ga^K - log(l + e^^)] - E log(l + e-^") 1 (6.1) 

I 6=1 j=l I 



where g = coxG, r = rankG, s = rank if. We introduce the notation GoH for the 
graph encoding of this TBA. This "product" is of course non-commutative, as one 
can not exchange the role of particles and magnons in general. The graph alone is 
not sufficient to encode the TBA: one still has to specify the form of the i/* . The rule 
is to encode masses proportional to the Perron-Frobenius of G, for a single node in 
H , while all the other are zero (or, when possible, left and right movers in the usual 
way). The diagonal TBA explored by Klassen and Melzer P| correspond to GoAi in 
this notation. For larger H one has to indicate on which node k the mass terms must 
be put, we do that by adding an index k to the whole GoH symbol. The magnonic 
TBA's studied in the present paper correspond to {Ai o H)^-, for all k E H . What 
we have proved in sect. 4 amounts to the statement that, considering general TBA's 
of the form ( |6.1| ), the case Go Ai admits sensible solutions only ioi G = A, D , E ,T 
and analogously AioH allows only for H = A,D, E, T. Unfortunately, we were not 
able to find a similar classification for the general GoH, in any case the set of G, H 



running on all ADET is already extremely rich. In |Tl[] the case of G = ADE and 
H = A only has been explored. We expect the other cases to hidden some beautiful 
surprise p9| . 



If, along the same lines of sect. 3, we search for the Y-system corresponding to 
the TBA (13), we get (here = e'-) 



y: [e + ^) y: [o-j)= n(i + nio)f- n(i + (6.2) 

This system shows a periodicity Y^(6 + inP) = Y^(9), with P = that encom- 
passes and generahzes both cases analyzed in this paper. 
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We would like to conclude by mentioning an intriguing observation. The Y- 
system ( |6.2D shows a curious duality: by exchanging e ^ —e we go from the 
Y-system of the G o H case to that of if o G. This clarifies why the two cases 
encompassed in the present paper had similar Y-system, related exactly by a sign 
flip in the e's. More generally one can think of some relation existing between the 
models described by TBA's dual to each other. We are at present not able to give 
any clear statement on this subject, simply we notice the amusing fact that the 
series of tadpole diagrams T^oAi considered by Klassen and Melzer to describe the 
perturbation of M2,2n+3 minimal models by their 0i3 field, goes under the described 
operation into the AioT^ case. Now, if we consider the {AioTn)i TBA, this describes 
the perturbation of M2„+i,2n+3 by its (pn field. What is surprising is that the two 
models M2,2n+3 and M2„+i,2n+3 are known to be related by level-rank duality |]30[| . 



If this important property of CFT has a relation to the duality described here or 



not, has to be explored further |]29 
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